Ground-based detection of a vibration-rotation line of HD in Orion 
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Abstract 

The V = 1-0 R(5) line of HD at 2.46 /iin has been detected at the position 
of brightest hne emission of shocked H2 in the Orion Molecular Cloud. The 
flux in this HD line, when compared to that of the previously detected HD 
0-0 R(5) line at 19.43 /xm, suggests that, like the v=l levels of H2, the 
v=l levels of HD are populated in LTE, despite their much higher rates of 
spontaneous emission compared to H2. The higher than expected population 
of vibrationally excited HD may be due to chemical coupling of HD to H2 via 
the reactive collisions HD + H ^ H2 + D in the shocked gas. The deuterium 
abundance implied by the strengths of these lines relative to those of H2 is 
(5.1±1.9) X 10-^ 
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1. Introduction 

Until recently, observations of the deuterium abundance in the interstel- 
lar medium have been restricted almost entirely to ultraviolet and millimeter 
wavebands. Both wavebands have limitations: UV observations probe only 
regions of low extinction; millimeter measurements are uncertain because of 
the inaccurately known degree of chemical fractionation, which favors for- 
mation of deuterated species at low temperatures. In addition, analyses of 
measurements in both bands are plagued by uncertainties in the corrections 
for line saturation (in the non-deuterated species). 

Recently, two pure rotational lines of HD were detected by the spec- 
trometers on the Infrared Space Observatory (ISO), toward the Orion star 
forming region OMC-1, offering new opportunities for determining [D]/[H]. 
Wright et al. (1999) used the Long Wavelength Spectrometer (LWS) to de- 
tect the pure rotational R(0) J— 1-0 transition at 112 /xm toward the Orion 
Bar. Bertoldi et al. (1999) detected the 19.43 yum pure rotational R(5) {v—0, 
J— 6-5) line at H2 Peak 1, the brightest position of shocked hne emission in 
the OMC-1 outflow, with the Short Wavelength Spectrometer (SWS). Both 
measurements have led to new estimates of [D]/[H] in Orion, Wright ct al. 
finding (1.0 ±0.3) x 10~^ and Bertoldi et al. obtaining (7.6 ± 2.9) x 10"^ 

Both of these estimates have large uncertainties, stemming only in part 
from the low signal-to-noise ratios of the line detections. One source of uncer- 
tainty is the value of [HD]/[H] in the warm partially dissociated gas produced 
by C-shocks, where the endothermic reaction, HD -|- H -|- 418K — > H2 -|- D, 
selectively depletes HD. Bertoldi et al. (1999) estimate that the depletion 
factor at Peak 1 is 1.67 for gas in the range of temperatures where most of 
the HD V — R(5) hne emission occurs. The second is the lack of informa- 
tion about the populations of other energy levels. This is a serious problem 
at Peak 1, because (1) the post-shock gas has a wide range of temperatures, 
and (2) HD radiatively relaxes much more rapidly than H2, so that its v=0 
J— 6 level, which is 2,636 K above ground, might not be in LTE, unlike H2 
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lines of similar excitation. 

The value of [D]/[H] given by Bertoldi et al. (1999) includes both the cor- 
rection for depletion and an adjustment factor of 1.5, their estimate of the 
factor by which the population of v—0, J— 6 falls below that for LTE. This 
factor, which is highly uncertain, could be tested if other transitions of HD 
were detected. One possible source of additional information is the funda- 
mental vibration-rotation band of HD, which should be excited into emission 
by collisions in the post-shock gas. The fundamental band occurs at the long 
wavelength edge of the 2 fj,m window and short wavelength edge of the 3 /im 
window, where telluric absorption, mostly by H2O, is severe. However, the 
v—1-0 R(5) hue of HD at 2.459 yum is accessible from a dry, high altitude site 
such as Mauna Kea. The upper level of this hue is 7,747 K above ground 
(Herzberg 1950; Essenwanger & Gush 1984), and hence the intensity ratio 
of it to the 0-0 R(5) line observed by ISO can provide valuable information 
concerning the excitation of HD. The 1-0 Q(5) line of H2, with upper level 
energy 8365 K (Dabrowski & Herzberg 1984), is nearby at 2.455 /xm, so its 
strength can be easily compared to the HD line. 

2. Observations and Results 

A search for the 1-0 R(5) hue of HD was conducted on UT 1999 January 
20 at the United Kingdom Infrared Telescope (UKIRT), using the 1-5 iim 
spectrometer CGS4 (Mountain et al. 1990). CGS4's echelle was used with a 
slit of width 0.82 arcsec oriented east-west across OMC-1 Peak 1 for a total 
of 40 minutes of exposure time; an equal amount of time was spent on blank 
sky 5 arcmin east. The resolution was 16 km s~^. The spectrum of the star 
HR 2007, measured with the same instrumental configuration, was used to 
flux- and wavelength- calibrate (using telluric absorption lines) the spectrum 
at Peak 1. 

The resulting spectrum of a 0.82x7.28 arcsec^ strip at Peak 1 is shown 
in Fig. 1. The spectrum is dominated by the 1-0 Q(5) fine of H2, but the 1-0 



3 



R(5) line of HD is clearly detected at about a thousandth the strength of the 
H2 line, roughly as predicted by Hartquist et al. (1982). The fluxes detected 
in the above aperture were (1.06 ± 0.21) x 10"^^ W and (1.2 ± 0.3) x 
10~^^W m~^, respectively; uncertainty in the flux cahbration dominates the 
uncertainty in the H2 line flux and is about half the uncertainty in the HD 
line flux. The HD line is at the predicted wavelength to within 0.0001 /xm. 
Its intensity distribution, which is similar to that of the Q(5) line along the 
slit, solidifies the identification. The profiles of the two lines do not appear 
identical, but are the same to within the noise. 

3. Analysis 

3.1 Beam dilution 

To determine the excitation parameters of the HD at Peak 1, it is neces- 
sary to compare the two line fluxes. The 0-0 R(5) line was measured by the 
SWS in an aperture of 380 arcsec^, centered on Peak 1, whereas the 1-0 R(5) 
line was observed at UKIRT through a 6 arcsec^ aperture on the brightest 
part of Peak 1. Prom the H2 Q(5) fine, which was observed both by the 
SWS in a large aperture (Rosenthal, Bertoldi, & Drapatz 2000) and UKIRT 
in the small aperture, we conclude that the observed HD 1-0 R(5) surface 
brightness must be decreased by a factor of (2.1±0.4) to properly compare it 
with the average surface brightness of the pure rotational HD line measured 
in the larger aperture. 

3.2 Column densities 

The surface brightness of an optically thin vibration-rotation line can be 
converted into a column density via N(v,J) = (47r/hc) (IA/A)10°-^^^), where 
I is the observed line intensity, A is the Einstein coefficient, and A^ is the 
extinction, assumed to be 0.78 mag at 2.455 ^m. (see Bertoldi et al. 1999, 
Rosenthal et al. 2000). Bertoldi et al. found an HD column density in the 
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v^O, J=6 level, N(0,6), of (3.0±1.1) x 10^^ cm-^ at Peak 1. After correcting 
for the difference in surface brightnesses as described above, the 1-0 R(5) hne 
yields a column density N(l,6) = (5.1±1.3) x 10^^ cm~^ averaged over the 
large beam. This is a factor of 60 lower than the column density in the 
v=0, J=6 level, demonstrating that most of the HD is not vibrationally 
excited. However, as shown in Fig. 2, the two populations match the relative 
populations of H2 levels (Rosenthal et al. 2000; Bertoldi et al. 1999) of 
similar energies. This is somewhat surprising, as one might have expected 
that the higher HD levels would have a sub-thermal excitation, due to their 
much higher radiative transition rates (60 times higher when comparing the 
rate for HD(1,6), 5 x 10""^ s~^, with that for H2(l,4), which has the same 
excitation energy). Thus one might have predicted that the population of 
the higher energy level of HD would fall far below the scaled H2 curve. We 
return to the discussion of HD excitation below. 

4. Deuterium abundance 

Assuming that the lower energy states of HD are populated in LTE (i.e., 
in the same way as the levels of H2, as suggested by Fig. 2), and after correct- 
ing for the depletion of HD, we find [D]/[H] — (5.1±1.9) xlO~^, where the 
uncertainty does not include the uncertainty in the depletion. This abun- 
dance ratio is 1.5 times lower than that given by Bertoldi et al. (1999), 
who had corrected for their assumed non-LTE level population by that fac- 
tor in addition to correcting for depletion. The abundance is the lowest yet 
determined for deuterium. However, within the uncertainties the value is 
consistent with the recent value of 7.4l};3 x 10~^ derived by Jenkins et al. 

(1999) toward 6 Orionis. We note that Jenkins et al. and Sonneborn et al. 

(2000) find that the deuterium abundance varies significantly along different 
lines of sight. 
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5. Excitation of HD 

The excitation of the HD (0,6) and (1,6) levels appears to mimic that 
of similarly excited levels of H2. If this is generally true for HD, then there 
must be a way of maintaining HD level populations in spite of their radiative 
transition rates, which are substantially higher than those of H2. While it is 
believed that H2 levels are in LTE as a result of standard coUisional relaxation 
with H2 and H, post-shock densities would need to be considerably higher 
than 10^ cm~^ in order for simple inelastic collisions to maintain HD in or 
near LTE. 

We have already noted the exchange reaction, HD + H ^ H2 + D, which 
tends to deplete HD in the warm, partially dissociated post-shock gas of 
OMC-1. Our analysis of this reaction, which has been studied previously 
by Rozenshtein et al. (1985), Zhang & Miller (1989), Gray & Balint-Kurti 
(1998), and Timmermann (1996), indicates that at Peak 1 the fraction of 
vibrationally excited HD, HD* , is significantly enhanced by this reaction. De- 
tails of our analysis are provided in Ramsay Howat et al. (2002). To sum- 
marize, we find the following. 

• H2(v=0) -|- D — > HD* -|- H is the dominant channel for "chemically" ex- 
citing HD vibrationally, with a rate coefficient of 1.24 x 10~^^ cm~^ s~^. 
Vibrationally excited H2 has an order of magnitude higher rate coef- 
ficient for production of HD*, but only a very small fraction of H2 is 
vibrationally excited (see below). 

• Formation of HD* is faster than its relaxation to the ground vibrational 
state when n(H2) > 4 x 10^[n(HD)/n(D)][ra(HD*)/n(HD)] cm-^. In 
this case collisions of HD* with H2 and other species will lead to an LTE 

population of HD*. The deuterium fraction, 'ri(D)/n(HD) is determined 
by detailed balance of the HD-H2 exchange reaction, and the above 
equation simplifies to n(H) > 2 x 10^ [n(HD*)/n(HD)] e*^^/^^^K 

• The ratio n(HD*)/n(HD), which is needed to evaluate the above in- 
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equality, varies greatly over the range of post-shock temperatures where 
hne emission occurs. For a reasonable assumption of 0.01 for the 
mean value of this ratio, the formation of HD* by reactive coUisions 
is faster than its radiative decay to the ground vibrational state when 
n(H) > 2 X 10^ cm-3. 

In the post-shock gas n(H) is believed to be well in excess of 10^ cm~^. 
Thus, we tentatively conclude that at Peak 1 partially dissociative shocks 
strongly couple the excitation of HD to H2. Because of the much larger 
abundance of H2, the HD effectively becomes part of the H2 level system. 

To further test the importance of this excitation mechanism, detections 
of additional lines of HD are needed, in order to determine if the populations 
of their upper levels also fall on the scaled H2 curve in Fig. 2. 
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Figure Captions 



Figure 1 - Spectrum of a 0.82" x7.28" (NS x EW) area of at the center 
of OMC-1 Peak 1 near 2.46 /im. The H2 and HD hues are indicated. The 
assumed continuum, used to estimate the flux in the HD hne, is shown by 
the dashed hne. 

Figure 2 - HD excitation diagram, with column densities plotted versus 
energy of upper level. Triangles are pure rotational transitions; diamonds 
are vibration-rotation transitions. Most measurements are upper limits. The 
curve is the fit to the H2 excitation diagram, scaled by a factor of 6.2 x 10~^ 
to pass through the HD 0-0 R(5) point. 
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